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lattice. In this paper a numerical approach based on the GAS algorithm is used to 
approximately enumerate states in the partition function of this model. The data are 
used to approximate the free energy in the model, from which estimates of the location 
of the critical point and crossover exponents are made. The critical point is found to 
be located at 

f 1.779 ± 0.003, in the square lattice: 
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1.306 ± 0.007, in the cubic lattice. 

These results are then used to estimate the crossover exponent (p associated with the 
adsorption transition, giving 


I 0.496 ±0.009, in two dimensions; 
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0.505 ± 0.006, in three dimensions. 

In addition, the scaling of these thermodynamic quantities is examined using the 
numerical data, including the scaling of metric quantities, and the partition and 
generating functions. In all cases results and numerical values of exponents were 
obtained which are consistent with estimates in the literature. 
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The adsorption of a linear polymer on an attractive surface is a conformational 
rearrangement of the polymer to a state where it explores conformations which remain 
near or on the surface. This is the so-called polymer adsorption transition, and the many 
models of this phenomenon (see, for example, reference |22|) remain a rich source of 
mathematical and numerical studies. Polymer adsorption is a phase transition [81, and 
the properties of the adsorbed polymer have been examined both experimentally (see, 
for example, references [7,13,32|) and theoretically 112]. These models include directed 

as well as self-avoiding walk models 1181, 


path models of adsorbing polymers [35 40 


and these have received considerable attention in the literature [ 191, using both rigorous 
methods |38| and numerical methods (for example, the Monte Carlo simulation of 
adsorbing self-avoiding walks |23|). 

In this study, a new Monte Carlo method for sampling adsorbing self-avoiding walks 
is proposed and implemented. In particular, the GAS algorithm |24| is generalised to 
sample adsorbing walks in the microcanonical ensemble, and the data obtained are 
analysed to estimate the locations of critical points, and the values of critical exponents 
and scaling of adsorbing walks, in the square and cubic lattices. The algorithm is 
related to the Rosenbluth algorithm |37| and to the GARM algoritm |36|. The GARM 
algorithm is related to the PERM algorithm 115,20], and flat histogram implementations 
of PERM |33| have been used to sample states from a flat histogram over state space 
in a variety of different models of interacting walks, including collapsing walks |34| and 
adsorbing walks |28|. 



Figure 1: An adsorbing self-avoiding walk in the half square lattice L+. The walk steps from the origin 
in the lattice, has length 46, and makes 6 visits to the adsorbing line dh+ (the boundary of L^). The 
weight of this walk in the partition function Z 46 (a) is a 6 . 


The GAS algorithm was introduced in reference |24| and was used in several studies 
as an algorithm to approximately enumerate walks or polygons |25,26]. However, it was 
not clear how to generalise the algorithm to sample states in models of interacting walks. 
In this paper our purpose is (1) to generalise the algorithm to a model of adsorbing walks, 
(2) to examine the behaviour of the algorithm by computing critical points and exponents 
of the walk, and to compare this to results found elsewhere, and (3) to use our data to 
examine scaling in adsorbing walks by computing critical exponents and examining the 
scaling of thermodynamic functions. 
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1.1. Adsorbing walks 

An adsorbing self-avoiding walk in the square lattice is illustrated in figure [l] Let h d 
denote the c/-dimensional hypercubic lattice and denote a unit length edge with endpoints 
x and y in h d by (x~y). The hypercubic half-lattice L+ is defined by 

L+ = {(x~y) g L d | x(d) > 0 and y(d) > 0}, (3) 

where x(d) is the d- th Cartesian component of x, and y(d) is the d- th Cartesian 
component of y. The boundary of is given by 

dh d + = {{x~y) e h d + | x(d) = 0 and y(d) = 0}, (4) 

and it is isomorphic to h d ~ 1 if d > 2. Please note that these definitions, and the 
definitions of additional functions and quantities, are listed in table [l] 

The number of self-avoiding walks of length n from the origin in L d is denoted by 
c n . The growth constant of the self-avoiding walk 116,171 is defined by the limit 

lim c n 1/n = p, d , (5) 

n—Foo 

and K d — \ogp d is the connective constant of the self-avoiding walk. This shows that 
c n — p n d + °^ n \ Self-avoiding walks from the origin in L+ are positive walks. The number 
of positive walks of length n is denoted by c+, and it is known that c+ = d n d + ° < ' n ' > |18|. 

If positive walks are counted with respect to the number of vertices in (these 
are visits ), then the walks are adsorbing walks. For example, the walk in figure [l] is an 
adsorbing walk with 6 visits. 

Let c+(v) be the number of adsorbing walks of length n from the origin in L+, 
with v visits to 9L^. The canonical partition function of adsorbing walks is obtained by 
introducing an activity a conjugate to the number of visits: 

n 

Z n( a ) = J2 C n^ 3V - ( 6 ) 

v=0 

When a is large, then Z n {a) is dominated by walks with a large number of visits, and if a 
is small (but positive), then Z n (a) is dominated by walks with a small number of visits. 

The finite size free energy of these models is computed from the partition function 
Z /7 (a) (see equation ([6])), and is given by 

fF n (a) = jj log ^ c+(v) a v . (7) 


The limiting free energy of the model is given by the thermodynamic limit in the model: 


T(a) 


lim 

n—>oo 


1 

n 


| °g^ c n + (' / ) 

v/ 


( 8 ) 


This limit exists (see reference 1181, and also, for example, reference |22|), and it is a 
convex function of log a with a singular point at a = a+ (which is the adsorption critical 
point in the model). For a < a+ the model is in a desorbed state, and for a > a+ the 


model is in an adsorbed state. It is known that at > 1 


and at < 


Me/ 

Md—1 


181, and 


Ha) 


= log p.d, 
> log p, d , 


if a < a+; 
if a > a+. 


(9) 
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Since jF(a) is a convex function of log a, it is differentiable for almost all a > 0, and it 
follows that the density of visits to the adsorbing plane is a-j^T(a) = 0 if a < a+ (this 
is the desorbed phase), and a^JF(a) > 0 if a > a+ almost surely (whenever jF(a) is 
differentiable). This is the adsorbed phase. In the desorbed phase the walk tends to 
make few returns to the adsorbing plane (walks of length n will return, on average, o(n) 
times to the adsorbing plane in the desorbed phase). Thus, a desorbed walk will tend 
to move away from the boundary into the bulk of L+. An adsorbed walk, on the other 
hand, is expected to have a positive density of returns to the adsorbing plane. This 
implies that the walk will remain near the adsorbing plane, and so have the properties 
of a walk which is stretched out in (d — 1) dimensions near 3L+ (and compressed in the 
d-th dimension). Separating these two regimes is the adsorption critical point a+. 

The finite size scaling of the free energy fF n (a) is given by 

T n (a) ~ log p, d + (a-a+) 2_a f(//(a-a+)), (10) 


where f is a scaling function, a is the specific heat exponent, and <p is the finite size 
crossover exponent. The exponents a and (p are related by the hyperscaling relation 

2 -« = ?- (n) 

The bulk entropy contribution to fF n {a) is \ogg d in L+ (where p, d is the growth 
constant of self-avoiding walks in the square lattice). Slightly redefining the scaling 
function, it is found that 

Fn(a) ~ log+ ^g(n 0 (a-a+)). (12) 

By plotting n(IF n (a)— log p, d ) against n^(a —a+), the function g can be uncovered (for 
n large and |a — a+| small). 

Taking derivatives of jF n (a) to log a gives the energy (density) £ n (a) and specific 
heat C n (a) of the model. The scaling of these quantities follows directly from equation 
©: 

£ n (a) ~ // _1 h e (n 4> (a — a+)), and C n (a) ~ n a<t> h c (n cp (a — a^)), (13) 


for some scaling functions h e and h c . In the limit as n —> oo, £ n (a) —>■ £(a) (the limiting 
energy density) and C n (a ) —>• C(a) (the limiting specific heat). Existence of these limits 
(almost everywhere) is a consequence of the convexity properties of the limiting free 
energy (see for example reference |22|). Physically, £(a) is the density of visits per unit 
length, and C(a) is the rate of change in £(a) as a function of changes in log a (it has a 
maximum at a+). 

For adsorbing walks it is thought that (p — \ in all dimensions d >2 |2,9|, and 


numerical evidence supporting this in dimensions lower than d — 4 (the upper critical 
dimension) are available in references |4 , 23 , 27 , 29 31 ] . If 0 = |, then a = 0, 
so, for example, the specific heat has scaling C n (a — h c (n (t> (a — a^)), and plotting 
measurements of C n (a) against the rescaled variable r = n^(a —a+) for small values 
of r should collapse the curves to a limiting curve (with some finite size corrections to 
scaling), exposing the scaling function h c . 
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The partition function has a more complicated scaling law (see, for example, 
equation (23) in reference |23|, or section 4.2.2 in reference [22]). In the high 
temperature (or small a) regime, the partition function scales as 

Z n (a ) ~ B A n 7t-1 /7x(^|t|) ^ , if a < a+, (14) 

where t — (a — a+) and A denotes the high temperature regime, and where h\(x) ~ 
| x |(7i-7t)/<#> anc j log ~ |a — a+| -1 ^ log p d . Putting a = 1, for example, and adsorbing 
constants and functions of t into B Xl give 

Z n (l) ~ B A n 7l_ Vd. since IF(a) = log p, d if a < a+, (15) 

and the exponent 'y 1 is the entropic exponent of half-space walks, namely c+ - n^p.% 
At the critical adsorption point a+, the above scaling is modified to 

Z n (a+) ~ B c n 7t “Vd. ifa = a+, (16) 

where 7 t is the entropic exponent associated with adsorbing walks at the critical 
adsorption point. The ensemble of half-space walks at the critical point has associated 
surface entropic exponent 7 S , and this is related to 7 t by 7 t = 7 S (see for example section 
9.1.3 in reference (22]). 

The scaling in the adsorbed phase is similar to the above, but now with different 
exponents 

Z n (a ) ~ B ro /7' Tt_1 /7 To (//|t|) k,'I w1/ \ if a > a+, (17) 

where h To (x ) ~ |x| (m+— anc | where the subscript r 0 denotes the low temperature 
(and large a) scaling. Since ^logZ„(a) = T'(a) (l + o(l)), it follows that logK + ~ 
[a-a+l-^^a). This, in particular, gives the scaling 

Z- s (a)-|a-a+| 1/ ^ (18) 

for the singular part of the free energy in the adsorbed phase (consistent with the 
hyperscaling relation (fTIj)). The scaling of Z n {a ) simplifies here to 

Z, 7 (a) ~ B To r7 7+_1 e n:F ^ a \ if a > a+. (19) 

The exponent 7 + should be that of adsorbed walks, and so given by 7 + = 7 ^ d_1 \ the 
entropic exponent of walks in one dimension lower. 


1.2. Organisation of the manuscript 


This paper is a report on two aspects of the Monte Carlo simulation of adsorbing walks. 
The first is the generalisation of the GAS algorithm to a model of interacting walks, 
and in particular, an implementation of this algorithm to achieve flat histogram sampling 
over state space of adsorbing square and cubic lattice walks. The second aspect of the 
paper is a report on the properties of adsorbing self-avoiding walks in the square and 
cubic lattices. The aim here is to verify the results obtained in the Multiple Markov 


Chain Monte Carlo study in reference |23|, and also to use the data generated here to 
test the scaling of the thermodynamic and metric quantities of adsorbing walks. 


The model of adsorbing walks is defined in section [lTT| and its partition function and 
free energy are discussed. The limiting free energy of this model exists, and its properties 
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Table 1: Short list of definitions 


Function 

h d 

14 

dh d + 

dd 

CnM 

Z n (a) 

H(a) 

Ha) 


Definition 

The hypercubic lattice 

The half-hypercubic lattice 

The boundary of (it is isomorphic is L^ 1 ) 

The growth constant of self-avoiding walks in L d 
The number of positive walks of length n and v visits from 0 in L+ 
Partition function of positive adsorbing walks of length n and activity a 
Finite size free energy: F n (a ) = ^ log Z n (a ) 

Limiting free energy: -T(a) = lim fF n {a) 


n—rOO 


S n (a), £(a) 
C n (a), C(a ) 
H n 

P+(6) 

at 

(p 

a 

7i 

u 

G(a, t) 
G/v(a, t) 
t+(a) 


Energy density and limiting energy density 
Specific heat and limiting specific heat 

The mean height of the endpoint of the walk (a function of a) 

The mean square radius of gyration of the walk (a function of a) 

The microcanonical density function (see equation (59)) 

The adsorption critical point 
The crossover exponent 
The specific heat exponent 

The half-space entropic exponent, see also 7 S (the surface exponent) 
The metric exponent 

Generating function of Z n (a) (see equation (66)) 

Truncated generating function (see equation (67)) 

The radius of convergence of G(a, t) 


have been examined elsewhere 118]; see, for example, reference |22|. The basic scaling 
relations for the free energy, energy, and specific heat were introduced above, and the 
scaling of the partition function was briefly reviewed. 

In section [2] the GAS algorithm |24| is reviewed and then generalised to interacting 
models. The algorithm normally has only one set of parameters (associated with the size 
of the walks), but it is shown here that introducing a second set of parameters (associated 
with the energy of the walks) can give an algorithm which samples effectively in both 
length and energy. It is shown that the algorithm can be tuned to give flat histogram 
sampling in the spirit of the PERM algorithm 115,201 (but without using enrichment or 
pruning of states). 

Numerical results for adsorbing walks are analysed in section [3] The location of the 
adsorption critical point is determined from the mean energy of adsorbing walks, giving 

f 1.779 ± 0.003, in the square lattice; 
at = < ( 20 ) 

I 1.306 ± 0.007, in the cubic lattice. 
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<■ 




Figure 2: Endpoint elementary moves. Appending an edge at the endpoint of the walk on the left is a 
positive elementary move. Removing the last edge on the right is the reverse of a positive elementary 
move. This is a negative elementary move. 


These results are then used to estimate the crossover exponent (p associated with the 
adsorption transition, giving 

{ 0.496 ± 0.009, in two dimensions; 

( 21 ) 

0.505 ± 0.006, in three dimensions. 

The microcanonical density function of adsorbing walks is determined as well, and shown 
to have properties consistent with the location of the critical points above. In addition, 
the specific heat of the model is determined, and it is found that it has scaling behaviour 
consistent with the estimates of the critical points above. 

It is also shown that the adsorption transition is seen in a change in the metric 
scaling of walks at the critical point. The desorbed phase is a phase of positive walks 
with the scaling properties of self-avoiding walks in a good solvent, while the adsorbed 
walk has the metric properties of a walk in one dimension less; that is, of walks adhering 
to the adsorbing surface. 

Scaling of the generating and partition functions are found to be consistent with the 
exact values of critical exponents determined elsewhere in two dimensions |2,9], and with 
the value (p — | in three dimensions 119,23,31]. Similary, the data are also consistent 
with estimates for the surface exponent 7 S ; this is seen particularly in the scaling of the 
partition function (see equation (fl6|)). 

The paper is concluded with a few brief remarks and a summary in section [4j 

2. GAS Sampling of self-avoiding walks 

The GAS algorithm is a generalisation of kinetic growth algorithms. It is designed to 
sample along weighted sequences in state space in such a way that the ratios of average 
weights of sequences ending in walks of length n and m are estimates of the ratio of 
the numbers of walks of lengths n and m. In this section I show how to generalise this 
algorithm so that it can be used to estimate the number of walks of length n and energy 
m. That is, the algorithm will be used to sample walks in the microcanonical ensemble. 

Let w — (tu 0 , u lt ..., cu n ) be a self-avoiding walk of length n from its source vertex 
cu 0 — 0 at the origin, to its terminal vertex u n , giving n steps {ujj-i, uij) for j — 1,2,, n. 
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Figure 3: A endpoint elementary move may increase or decrease the number of visits in a walk. By 
appending an edge at the endpoint of the walk on the left a new visit is created. Removing the last edge 
on the right is the reverse of this move, and it decreases the number of visits by one. 


The walk w may be made longer by adding a step (u n ,u n+1 ) to u n , or it may be made 
shorter by removing its last step. These two operations compose an end-point elementary 
move for sampling self-avoiding walks, as illustrated in figure [2} A positive elementary 
move is the addition of an edge to the endpoint of a growing walk. The reverse of a 
positive elementary move is a negative elementary move, namely the deletion of the least 
edge in a walk. Notice that every positive move is immediately reversible by a negative 
elementary move. 

Endpoint elementary moves have been used widely in the simulation of self-avoiding 
walks (for example the Rosenbluth algorithm |37|, and the Beretti-Sokal algorithm |6|). 
In what follows the discussion will be restricted to endpoint elementary moves; however, 
the algorithms generalise directly if other elementary moves, such as BFACF elementary 
moves 


or generalised atmospheric moves [361, are used instead. 

As an example, consider the elementary move in figure[2j which is an implementation 
of an endpoint elementary move on a self-avoiding walk in the positive half-lattice IUjr 
(where w(2) is the y-coordinate of w). The set of lattice edges in which may be 
appended to the walk \n to extend it by one step is the positive atmosphere of w, and 
the number of edges in the positive atmosphere is denoted by a + (w). For example, for 
the walk on the left in figure [2| a + (w ) = 2. 

Similarly, the set of edges which may be removed from the endpoint of a walk 
to decrease its length by one, is the negative atmosphere of the walk. For endpoint 
elementary moves, the last step is always the sole negative atmospheric edge, so that 
the size of the negative atmosphere for endpoint elementary moves is always a~(w ) = 1, 
if w is not the trivial walk of length 0. 

An elementary move may change the energy of a walk. For example, in a model of 
adsorbing walks in iJf, the energy is the number of visits of the walk to the adsorbing 
line dlJf (the boundary of the half-lattice Lj). The positive move in figure [2] does not 
change the number of visits, but the move in figure [3] increases the number of visits (and 
so the energy) by 1. The negative move in figure[3]similarly reduces the number of visits 
by 1. 

A similar situation arises if a model of collapsing walks with energy given by nearest 
neighbour contacts between vertices in the walk which are adjacent in L, but not in the 
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Figure 4: A endpoint elementary move may increase or decrease the number of contacts in a walk. By 
appending an edge at the endpoint of the walk on the left two new contacts are created. Removing the 
last edge on the right is the reverse of this move, and it decreases the number of contacts by two. 


walk. This is illustrated in figure]?] the positive elementary move creates 2 new contacts 
in the walk, and so it changes the energy of the walk by 2 . 

More general elementary moves (for example, the BFACF elementary moves) may 
contribute to the atmospheric statistics a + and a~ in various ways, and may even give 
rise to neutral atmospheres which do not change the length of the walk (but which may 
change the energy of the walk). 

Thus, in what follows, let a+(w) be the size of the positive atmosphere of a walk 
w of length l n — \w\, and which changes the energy of w by v units. For example, the 
walk on the left in figure [ 2 ] has 3 q(m/) = 1 and a^(w) — 1 , and if this was a walk in 
the three dimensional half-lattice L 3 , then a^(w) — 3 and a^f(w) — 1. The walk on 
the left in figure]?] (with energy given by the number of contacts), has, in a similar way, 
a o( w ) — 2 , 3^(w) — 0 and aj( 1 / 1 /) = 1 . 

In exactly the same way one may define the size of the negative atmosphere of a 
walk w of length i which changes the energy by v, denoted by a~(w). 

The neutral atmosphere (an elementary move which does not change the length 
of the walk) of a walk w is similarly given by a°(w), if it changes the energy of a walk 
by v. If the endpoint elementary move in figure [ 2 ] is used, then a°(w) — 0 by default 
(since there are no neutral elementary moves implemented), but in general assume that 
a more general set of elementary moves (for example BFACF elementary moves) is used 
to sample walks, and in that case the neutral atmosphere may have positive size. 

2.1. Implementation of GAS-sampling 

Suppose that an elementary move is implementated on the state space of self-avoiding 
walks from the origin, and assume the implementation is irreducible (that is, the 
elementary move gives a connected graph on the state space of walks). 

Suppose that the sequence 

4>n = (w 0 , wi. w 2 _ w n _ w N ) (22) 

is realised after N steps and that the atmospheres of the states w n have sizes (or 
statistics ) a+(i/V„), a°(w n ) and a~(w n ). These elementary moves may be classified as 
follows, depending on whether they increase or decrease the lengths of the walks (or are 
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neutral), or whether they increase or decrease, or leave unchanged, the energy of the 
walk. This is done by defining atmospheric statistics as follows: Let the states w n be 


walks of length £ n 

— \ w n \. Define 





a n~ 

= XXM' 

° L n° ~- 

= ^3 v ( w n ), 

a n + 

= ^a v (w n ) 


y<0 


v=0 


v>0 

a°~ 

= ^ a °(^)' 

«° n ° = 

-- 

ct 0+ = 

: ^a°v(w n )] 


v<0 


v=0 


v>0 

a + ~ 

^ n 



= ^4( W n). 

a ++ 

= ^at( w n) 


v<0 


v=0 


v>0 


For example, a~~ is the number of negative elementary moves which also decreases the 
energy of the walk, and a~ + is the number of negative elementary moves which also 
increases the energy of the walk. The rest of the a's are similarly defined. 

With these atmospheric statistics defined, a rule needs to be constructed in order 
to realise the sequence <p N in equation ( 22 ). 

The endpoint elementary moves in figures [3]and[4]have the property that no positive 
elementary move can decrease the energy, and no negative elementary move can increase 
the energy. Moreover, there are no neutral moves amongst the elementary moves. Thus, 
assume that 


a n + = 0, a+ =0, a° n =0, a°° = 0 and a° n + — 0. (23) 


The algorithm can be modified appropriately to account for such transitions in models 
where this is not the case. This assumption leaves the following atmospheric statistics: 
{a““, a“°, a+°, a+ + }. That is, distinguish between negative moves which decrease the 
energy, or negative moves which leave the energy unchanged, or positive moves which 
leave the energy unchanged, and positive moves which increase the energy. 

Introduce parameters {fy, u } to control positive elementary moves on walks of length 
i and energy u, and which leave the energy unchanged (that is, the elementary moves 
contributing to a+°). Similarly, introduce the parameters { 7 ^} to control positive 
elementary moves which increase the energy on walks of length t and energy u. 

The transition probabilities of positive elementary moves which leave the energy 
unchanged will be proportional to if the state has length i and energy u. For 
example, since i — 16 and u = 2 in the walk in figure |2j the transition probability of the 
positive move in that figure is proporsional to /3 i 6 , 2 . Similarly, the positive elementary 
move in figure [3] increases the energy; so here the transition probability is proportional 
1° 7 i 6 , 2 i instead. In figure [4] the positive elementary move also increases the energy, and 
so its transition probability is proportional to 7 i 6 , 4 - 

Thus, if w n is the current state (of length l n and energy u n ), and w n+1 is the next 
state (of length l r+ 1 and energy u n+ 1 ), then define the change in length by A n — i n +i~ln, 
and the change in energy by 5 n = u n+1 — u n . Notice that A„ = ±1 for endpoint elementary 
moves, and that 5 n = ± 1 , or 5 n — 0 , for the model of adsorbing walks in figures [ 2 ] and [3] 
(but these quantities may take on other values in the model in figure [4]). The transition 
probabilities are chosen such that 


if A n 
if A„ 


+ 1 and 5 n — 0; 
— 1 and 5 n — 0, 


P r {w n ->• w n+ 1 ) oc 


Plu 

1 


(24) 
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where i — i n and u — u n are functions of n. In the case that the energy is changed, then 
the transition probability is constructed such that 

if A n = +1 and 5 n > 0; 


P r (w n —> W n+ i) (X 


Itu, 


if A n = — 1 and 5 n < 0. 


(25) 


Normalising the transition probabilities gives 

f (3 i, u 


P r (w n —> W n+ i) — < 


Ot~~ + °T° + <Xt°Pl.u + OL++ 7^,, 

_ Itu _ 

a~~ + a ~° + ol^°Pc,u + ol+ + 7 ^, 
1 

ot~ 


if A n = 1 and 5 n — 0; 
if A n = 1 and 5 n > 0; 
if A n = —1 and 6 n < 0; 


(26) 


+ a-° + a+°P t , u + a++7(,i/ ’ 

The probability for any type of move can be explicitly computed for any given 
state w n (of length l and energy u ) by computing the a's. For example, the 
probability for a positive elementary move which leaves the energy unchanged is 
&~!i 0 Pe,u/{oc~- +a+°/3^ iU +a+ + 7£ iU ). For the state on the left in figure [2] this becomes 
016 , 2/(1 +/3i6 ,2 + 7i6, 2 ) since a““ = 1, a+ 6 ° = 1 and a++ = 1. 


2.2. GAS-weights 


Let w 0 be a starting state (possibly the walk consisting of a single vertex at the origin, of 
length — 0 and energy u 0 — 0). Implement the endpoint elementary moves on w 0 by 
computing its atmosphere recursively and updating it using the transition probabilities in 
equation (26), starting at n — 0. This generates a Markov Chain of states in a sequence 
( p N (see equation (22)). 

Assume that the parameters {P n ,u,ln,u} are known and fixed, so that the sampling 
can be implemented by simply computing the transtition probabilities and selecting 
positive and negative elementary moves with appropriate probabilities. 

The probability of the sequence <p N is given by 


N—l 


Wn) n 


n =0 n 


+ a F° + 0i+°Pi.u + o.+ + ^i u 


n' a- ii/"- 


(27) 


where the primed product is over all f-^Ptu for transitions through positive elementary 
moves leaving the energy unchanged, and the double primed product f| // is over all the 
where the transition is a positive elementary move increasing the energy. In this 
expression the L and u are functions of n, m, and k in each of the products. 

A weight W(<p N ) will be assigned to the sequence (p N . In order to compute the 
weight, define 


<r(jJ+ 1) = w j+ 1 ) 


— 1, if Aj = +1 and 5j — 0; 

+ 1, if A j — —1 and 5j — 0; 

0, otherwise. 


(28) 










Microcanonical Simulations of Adsorbing Self-Avoiding Walks 
Similarly, define 


12 


P(JJ+ !) = P( w j W/+i) = { +1. 

0, 


if Ay = +1 and <5y > 0; 
if Ay = — 1 and <5y < 0; 
otherwise. 


(29) 


That is, the function o(jj + l) tracks the negative and positive moves along <p N where 
the energy is not changed, and the function p(j,j + l) tracks the negative and positive 
transitions along <fi N where the energy is also changed. 

Assign the weight 


w(<p N ) = 


Oi 0 +g 0 °+Q!^ 0 /3o+a^ + 7o 

a N~ + a N°+ a ti°PN+al + 'y N 


N-l N—l 

n^iKT 11 

j =0 j =0 


(30) 


to the sequence <p N . 

The expected value of the weight over sequences of length N from state w 0 to state 
w N is 


(W(w 0 —» w N )) N — 


Cp-.W 0 ^W N 




(31) 


Inserting equations (|27]) and fl30| ) in this, and simplifying, gives 

N 1 

E l—r _L i —r n i— r-nn 


+a, °+ay°/3^ u +a, ++ 7f, u 


(32) 


{l/l/(w 0 —> w n )) n — I f 


(p'.WQ^Wi\l j= 1 J J J 

where the product is over the (3 tiU for transitions through negative elementary 
moves leaving the energy unchanged, and the product Y\ nn ' s over a H li.u where the 
transition is a negative elementary move decreasing the energy. As before, the i = i(j) 
and u = u(j) are functions of j (in other words, functions of the states Wj in the sequence 
(p). 

Reverse all the sequences in equation (32) so that the starting state is w N and 
the final state is w 0 . Under this reversal all negative elementary moves become positive 
elementary moves and vice versa. That is, equation (32) becomes 

N 

I -■—r n -m — w- nn 

• I Pi,, n ^ ( 33 ) 

~J 1 “J f'-'t.L/ 1 “j 1 iL,U 

where, as before, the product Yl' 1 is over a H the /3*, 0 for transitions through negative 
elementary moves along the reverse sequence ip leaving the energy unchanged, and the 
product FTis over a II li.u where the transition is a negative elementary move along the 
reverse sequence ip decreasing the energy. 

For example, consider the model of collapsing walks in figure [4] and suppose the 
sequence <p is realised, where 

*=(.-►_ “□). 

W 0 Wi W 2 W 3 W/\ W 5 W 5 

The probability of this sequence is 

p _ Po,Q Pl,0 _ J2X) _1_ J2S) _ fe,l 

r\V) 4/3o,o (l+3/3i,o) (1+2/32,o+ 72 ,o) (1+2/33,l) (1+2/32,o+ 72 ,o) (1+2/33 ,i) 


a +a u +ocf + 'Ye, 

ir.w N —>W q j=1 J 1 1 1 


( 34 ) 
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The weight of 0 can similarly be computed from equation (30). This gives 

W((p) = ! x (^ 0 , 0 /^i,o 72,o T^i 72,0^3, 1 ) ■ (35) 

The consequence is that 

P(0)kl/(0) (l+3/3 4 1 ) (1+2/3 3 i) (1+202,o+72,o) (1+203,i) (1+202,o+72,o) (l+30i,o)' (36) 

and this is the probability that a sequence ip, starting in the state w 6 and terminating in 
the state w 0 , is realised by the algorithm. 

The same observation is true generally for equation (33): The summand is the 
probability that a particular sequence ip from state w N to w 0 is realised by the algorithm, 
and the summation is over all such sequences ip. That is, {W(w Q -+ w N )) N — P r (w N —>■ 
i+o) is the probability that the algorithm realises a sequence ip of length N from state 
w N to state w 0 . 

The sequence ip is a Markov Chain, and if it is aperiodic and irreducible, and 
if w 0 is a recurrent state, then asymptotically (for large N ) the probability that the 
sequence hits the state w 0 is positive and independent of the starting state w N . 
That is, P r (w N —* i+ 0 ) —* C(w 0 ) > 0 as N —» oo where C(w 0 ) is dependent on 
the parameters of the algorithm, and independent of w N . Thus, the average weight 
(W(w 0 -+ w N )) N —> C(i+ 0 ) as N —> oo. Summing (W(w G —> w N )) N over all the states 
w N of length n = £(i+ w ) and energy u, shows that the average weights of sequences of 
length N ending in walks of length n and energy u is 

{Wn, u )i\i= z (I/1/(m/o + %))« + 1 c( u ) 6(wo), (37) 

w N : n,u 

where the summation is over all walks of length n and energy u (and c n (u ) is the number 
of walks of length n and energy u). 

Taking ratios of average weights give 


(W n , u ) 


N 


-+ 


c n (u) 


(38) 


(W m , v ) N c m (v)' 

That is, if c n (u ) is known for some values of n and u, then the ratios of average weights 
can be used to estimate c m (v). 


2.3. Sampling with GAS 

The algorithm is implemented by choosing a starting state w 0 and then sampling along 
a sequence <p of length N. The weight is updated along <p, and collected into bins for 
walks of length n and energy u. Once the sequence is completed, then the average of 
each bin is computed, giving the average weight W n u of walks of length n and energy u 
seen along <p. 

This process is repeated M times, so that M sequences of length N, denoted 
by (0i,02. ■ ■ ■, (Pm) are realised, and for each sequence 0 7 the average weight W$ u is 
calculated. The estimated average weight is computed over the M sequences: 

M 

i 


(39) 
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Table 2: Numerical estimates of @ n u for adsorbing walks in L+ 

n\u 012345678 9 10 

0 i.ooo 

1 0.333 1.000 

2 0.428 0.250 1.001 

3 0.368 0.499 0.200 1.000 

4 0.387 0.380 0.416 0.166 1.000 

5 0.373 0.396 0.427 0.429 0.142 1.000 

6 0.386 0.375 0.377 0.359 0.438 0.124 1.000 

7 0.377 0.392 0.384 0.401 0.348 0.445 0.111 1.000 

8 0.383 0.378 0.388 0.356 0.368 0.327 0.449 0.100 0.999 

9 0.378 0.386 0.385 0.395 0.370 0.368 0.308 0.454 0.0908 0.999 

10 0.381 0.379 0.383 0.373 0.368 0.352 0.364 0.289 0.458 0.083 0.999 


The estimated weight [W nu \ e N st M is an estimator of {W niU ) N , and so as M —> oo and 
N —f oo, it is expected that [W n u \ e N st M —> c n (u ) C(w 0 ). Taking ratios give 

inos. fW ,, m 

The sampling requires that both N and M are sufficiently large, and there is a trade¬ 
off between these quantities. M should be large enough to have sufficient independent 
estimates of the weights to compute sound averages, and N should be large enough to 
have sufficiently long sequences to have sampled large enough regions of state space. 

There remains the additional issue of the GAS parameters / 3 n , u and 7 „ iU . These can 
be estimated using training runs prior to the simulation. Best results are obtained when 
the sampling is flat, so that states of size and energy {n, u} are sampled uniformly in 
{n, u}. This is best achieved when the GAS sequences are random walks on n and u. 
Thus, the probability of a positive move should be, on average, equal to the probability 
of a negative move. A good choice for (3 n u is 



(41) 


where is the average negative atmosphere which does not decrease the energy, 

and (a+°) is the average positive atmosphere which does not increase the energy, of 
walks of length n and energy u. 
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Table 3: Numerical estimates of 7 n u for adsorbing walks in lp_ 



A similar argument shows that a good choice for 7 n u is 


7 n.u = 


a 


n,u 


a 


++ 


> 

> 


(42) 


where (oc~~) is the average negative atmosphere which decreases the energy, and («++) 
is the average positive atmosphere which increases the energy, of walks of length n and 
energy u. 

Computed values for p nu and 7 „ it , for a square lattice adsorbing walk model in figure 
[2] are given in tables [2] and |3] for 0 < n < 10 and 0 < u < 10. These data are chopped 
at an accuracy of three decimal places. 

Observe that p nill — 1 (when rounded) and 7 nin = 1. This is expected, since there 
are exactly two states if n > 0 and n — u (a completely adsorbed walk which never leaves 
the adsorbing line). 

With these values of P nu and 7 n u the sampling is reasonably flat, as shown in table 
[ 4 ] The data in table[ 4 ] is the number of times a sequence of length 10 6 visited states of 
length n and energy u. The maximum length was set at n = 500, so that the number 
of pairs ( n , u ) is 124750; thus, the expected number of visits to states of length n and 
energy u is roughly 8000. The data in table [4] are spread around this number, and 
the distribution is a reasonably flat histogram. There are exceptions for data along the 
diagonal where a larger number is seen. The explanation for this is that the sequence 
can only visit states along the diagonal if it starts in (n , u ) = (0,0) and then stay on 
the diagonal (doing a random walk on states of length and energy both equal to n ) - 
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Table 4: The frequency of states sampled to n = 10 and u — 10 for adsorbing walks in 


n\u 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

8658 











1 

8616 

13208 










2 

9642 

7711 

13261 









3 

8971 

9622 

7144 

13012 








4 

8919 

9576 

9261 

6528 

12883 







5 

8957 

8947 

9501 

9388 

5998 

13072 






6 

8982 

8736 

9220 

9019 

9281 

6055 

13263 





7 

9059 

8562 

9026 

9095 

8711 

9733 

5812 

13278 




8 

8996 

8398 

9030 

8844 

8568 

8914 

9493 

5728 

13362 



9 

9169 

8162 

8788 

8772 

8579 

8843 

8417 

9689 

5678 

13204 


10 

9134 

8182 

8580 

8434 

8559 

8898 

8390 

8221 

10101 

5477 

12883 


the elementary moves chosen in the simulation do not include neutral moves, and so the 
sequence cannot enter the diagonal except at n = 0. 

Observe that since c+(n) = 2 for n > 0 in this model, that these states are rare, but 
are sampled frequently by the algorithm. This is an example of rare event sampling, where 
an algorithm spends significant time sampling rare states in the tails of a distribution in 
order to get good estimates of microcanonical quantities. 

Simulations were performed by collecting data over 500 realised sequences, each of 
length 10 9 . Numerical estimates of c n<u are shown in table [5] The data in tableware 
rounded to the nearest integer. These data can be compared to exact counts to verify 
the algorithm and its implementation. 

Partition functions (see equation ([6])) can be directly estimated from the 
microcanonical data in table [5j for example, for n — 6 the partition function is 
approximated by 


Z 6 (a) = 131 + 106a + 56a 2 + 28a 3 + 14a 4 + 2a 5 + 2a 6 (43) 

where a = e v i kT is a Boltzmann factor (•v is the interaction energy associated with a 
single visit). 

Additional data were collected in the square and cubic lattices. For example, the 
average mean square radius of gyration of walks of length n and energy u were also 
determined, as was the average height of the endpoint of the walks. In addition, data 
were also collected on the average positive and negative atmospheric statistics from 
which estimates of the f3 n<u and ^ n u were made (see equations (|4lj) and (42)). 
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Table 5: Estimates of c+(u) to n = 10 and u = 10 for adsorbing walks in L+ 



3. Numerical Results 


3.1. Adsorbing walks in the square lattice 

The finite size free energy fF n {a) was determined from the data and is plotted in figure [5] 
for walks of lengths n — 50 N with = 1,2,3,..., 10. LF n {a) is a function of the combined 
variable r = n 1//2 (a — a+) (see equations (10) and (12), and note that the crossover 


exponent is 0 = |). Plotting the free energy tF n {a) against r should collapse the data to 
a single underlying curve (near the critical point a+; that is, for small values of r) which 
exposes the scaling function g in equation (12). This is displayed in figure [6] 

The energy £„(a) and specific heat C n (a ) (see equation (13)) can be determined by 
differentiation of T” n (a) and is plotted in figure [7] against log a and in figure [8] against r. 
These curves clearly show the adsorption transition in the model at a critical value of a. 

Closer inspection of the specific heat curves in figure [7] shows that they intersect 
each other close to a fixed point. To the left of this point the curves decrease with 
increasing n to zero, and to the right of this point the curves increase with increasing n. 
The common point of intersection is located approximatedly at the critical adsorption 
point a+ in the model. In figure [9] the specific heat curves are magnified in a region close 
to the point where they intersect. 


3.1.1. Location of the critical point a+; In order to determine a+, consider the finite 
size energy density £ n (a) = a^T'(a) per edge in the model. It is known that (see, for 
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Figure 5: The finite size free energy T n {a) of adsorbing positive walks as a function of a for adsorbing 
walks in Lj ■ The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to 
black when n = 500 in steps of 50. For small values of a tF n {a) converges to log ^ 2 , but for a large 
T n {a) > log H 2 - 


example, equation (13) and reference 1231) 
£ n (a) = n 0-1 h n (j) 


where h n is a scaling function. From this equation, construct the ratio 
log (nS n (a)) _ (p log n + log h n (r) 
log (m£ m (a)) (p log m + log h m (r) 


(44) 

(45) 



Figure 6: The scaled free energy n(tF n — log ^ 2 ) as a function of r = n 1 / 2 (a— ap), where ap = 1.78, 
for adsorbing walks in hp. The lengths of the walks varied from n = 50 (yellow) with colours increasing 
in hue to black when n = 500 in steps of 50. The data collapse for small |r| (this is the critical scaling 
regime), and will approach a limiting curve as n —> 00 , also for large r, given by the limiting free energy. 
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Figure 7: The energy £ n (a) and specific heat C n (a ) (see equation (131) as a function of log a for adsorbing 
walks in Lj. The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black 
when n = 500 in steps of 50. 


If a — a+, then r — 0, and the above simplifies to 
log(ng n (a+)) = 0 log n + log h n (0 ) 
log (m£ m (a+)) 0 log m + log h m (0 )' 

For large values of n and m, log h n (0) and log h m (0) approaches the same constant C, 
so that the above becomes 

log(/7 £„(a+)) 0 log n + C 


\oq(m£ m (a+)) 0 log m + C 


+ small correction. 


(47) 



Figure 8: The rescaled energy density (left) and the specific heat (right) plotted as a function of 
t = n 1//2 (a — a+) near the critical point a+ for adsorbing walks in Lj- The curves approach a limiting 
curve with increasing n. The lengths of the walks varied from n = 50 (yellow) with colours increasing in 
hue to black when n = 500 in steps of 50. The data collapse for small values of |r| (this is the critical 
scaling regime which contains the critical point). For larger values of |r| the curves will approach a limiting 
curve as n —> oo given by the limiting free energy. 
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Figure 9: A magnification of the point close to where the specific heat curves in figure [7] intersect, for 
adsorbing walks in Lj ■ The lengths of the walks varied from n = 50 (yellow) with colours increasing in 
hue to black when n = 500 in steps of 50. 


When C is small compared to log m, the right hand side may be expanded to obtain 


log (/i £ n (a+)) 
log (mS m (a+)) 


log n f C C log n 
log m + [4> log m log 2 m 


C 2 

<p 2 log 2 m 


+ small correction. 


(48) 


The signs of the two terms in square brackets are opposite, and for n and m not too far 
apart, these terms grow about at the same rate as the last term ^ m . That is, the 
approximation 


log(n£ n (a+)) 
log (m£ m (a+)) 


log n 
log m 


Co 

<t> 2 log 2 m 


+ small correction, 


(49) 


where C 0 is a constant, should be accurate at the critical point. Divide both sides by 
^ to obtain 

log m 


Pn,m(at) = 


log(ng„(a+)) log m 
log (m£ m (a+)) log n 


= 1 - 


c 2 

d 2 log n log m 


+ small correction, (50) 


Solving for the critical point by inverting P n:m gives the solution a+ m as an approximation 
of the critical point, namely = P~^( 1 — ^ | og C ,)| ogm Tsmall correction) which may be 

expanded to 

= Pn.M - log n^og m + smal1 correctlon ' ( 51 ) 

where C 1 is a constant. That is, for given values of n and m, an estimate of a+ can be 
obtained by determining the solution of P n , m (a ) = 1, or a+ m is estimated obtained by 
solving for a in 

log(n£ n (a)) log m 
log (m£ m (a)) log n 


1 . 


(52) 
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Figure 10: Extrapolating the intersections between specific heat curves for adsorbing walks in L.f. These 
data points are the intersections of the specific heat for n and n+ 100, for n even and n e [44,400]. The 
data (except for some data points at the largest values of ri) line up along a line if plotted against 1 /yT. 
The location of the critical point can be estimated by extrapolating the line to the left vertical axis; this 
gives the estimate in equation ((551). 


In the above, m was put equal to n— 100, and n was assigned values starting at n — 200 
to n — 500 in steps of 1. The estimates 3+ = a„ J1 _ 100 showed no systematic dependence 
on n, and a simple average over all n e [200, 500] gives the best value 

a+ = 1.7788 ± 0.0029. (53) 


The confidence interval is obtained by doubling the square root of the variance of 


the estimates 3+. This result compares well with the result in reference 


namely 


= 1.77564 (obtained by the exact enumeration of adsorbing walks), and also with 
3+ = 1.759 ± 0.018 in reference |23| (obtained by using a Multiple Markov Chain 
implementation of the Berretti-Sokal algorithm |6|). 

The estimate (53) can be used to determine the crossover exponent 0. By equation 
(13), the specific heat scales as C n (a+) ~ n a ^/i c (0) when a — a+. That is, an estimate of 
a<p is obtained by computing log ^ n g ^ c ~ • Computing this for 100 < n < 500, and taking the 
mean as the best estimate (and estimating a confidence interval by doubling the square 
root of the variance of the estimates), gives a0 = —0.0091 ± 0.0162. Determining 0 
by using equation (JTTj) then gives the best estimate for 0: 

0 = 0.4955 ± 0.0081. (54) 

This result is consistent with 0 = | |2,4,9j, and compares well with other estimates in 
the literature (for example, 0 = 0.501 ± 0.014 in reference |23|). 


3.1.2. The critical point 3+ and the specific heat C n (a): The best estimates above may 
be compard to estimates obtained from the specific heat curves in figure[7J These curves 
intersect each other near a+, and the region containing the intersections (in figure [9]) is 


magnified in figure 10 
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In general the location of the intersection between C n (a ) and C m (a) is a function 
of n and m. The location of the critical point a+ can be estimated by extrapolating 
this dependence. Consider for example the intersections between the curves C n (a) and 
C n +ioo(a)- The location of these intersections are plotted against ^ in figure 
n — 2 N and N e [23,200]. The data lie along a straight line, except for a 
at the largest values of n (where the data is more uncertain). The best line through the 
data can be extrapolated to its intersection with the vertical axis (where ^ = 0). This 
gives an estimate of a+ as being located in the interval [1.77,1.79]. Using a linear least 
squares model for n > 10 gives the (extrapolated) estimate a+ ~ 1.7839. By examining 
the spread of the data in figure[lOj a confidence interval can be estimated. The result is 


Toj where 
ew points 


at = 1.784 ±0.010. 


(55) 


This estimate is consistent with the best estimate in equation (53). 

An alternative approach to determining the critical point at, is to consider the scaling 


of the specific heat in equation (13). Taking ratios for n and m, and then logarithms, 
give 


lo 9 (fell) ~ l°g (m) T log 


/?c(n 0 (a-a+)) 

hcim^ia—at)) 


(56) 


Observe that the last term is zero when a — at. Since, in addition, a = 2—} = 0 in 
this model, this shows that an estimate of the critical point is given by the solution of 

'°S (S§) = 0. (57) 

Solving this for 150 < n < 500 (and n a multiple of 10), and for m — n—k where 
k e {10, 20 ,..., 100} gives a large collection of estimates of a+, with mean 


a+ = 1.762 ±0.016. 


(58) 


The confidence interval is one-half of difference between the maximum and minimum 
estimates of a+. This estimate is slightly less than,but still consistent with, the results 
in equation (53) and equation (55). Recall that it is also predicated on the assumption 
that a = 0 (or 0 = |). 


3.1.3. The microcanonical density function: The microcanonical density function of 
visits in adsorbing positive walks is determined from the microcanonical data in the 
model, and is given by 

P + (e) = lim (c+(|_enj)) 1/n = lim P„ + (e), (59) 

n—¥ oc n—>oo 

where P+(e) = (c+QefiJ)) 1 /' 7 is a finite size approximation to P + (e). Existence of the 
limit is known (see for example reference [22]), and log P + (e) is a concave function of 

e. 

P + (e) can be estimated by interpolating the finite size approximations P+(e) and 
then extrapolating to n = oo by fitting a least squares model to the data. In figure[TT] 
the data for the extrapolated function P + (e) is plotted together with P n + (e) for n — 100 
and n — 500. 
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Figure 11: The microcanonical density function P + (e) of visits of adsorbing square lattice walks. These 
data the finite size approximations P+(e) to P + (e) for n = 100 and n = 500, as well as the extrapolated 
estimate of P + (e). The right derivative of log P + (e) at e = 0 is an estimate of the location of the critical 
adsorption point a+ in the model (see, for example, reference 


22 


A least squares fit of a quadratic to log R + (e) for e e [0, 0.1] gives the log R + (e) ~ 
0.97007 —0.58190e — 0.13030e 2 , and by taking the right derivative and then taking 
e —y 0 + , an estimate for the critical point is obtained: 


a T - 

c ~ 


1.789. 


(60) 


This is close to the estimates obtained in equations (53) and (55), showing consistency 
in the data and the analysis above. 

The free energy jF(a) is the Legendre transform of log R + (e). This may be estimated 
by fitting a polynomial to log R + (e). If a cubic polynomial in e is fitted to log R + (e) for 
0 < e < 0.5, then the estimated free energy for a > at is approximately 

F(a) « 1.0376 - 0.1153 log a - (0.3306-0.5693 log a)V-1.3603+2.3422 log a. (61) 


The critical point can be estimated as that location where the square root in the above 
is zero. This gives a+ ~ 1.787. Similarly, the factor (0.3306 — 0.5693 log a) vanishes 
when at ~ 1.787. 


3.1.4. Metric data: The mean square radius of gyration R 2 , and the mean height H n 
of the endpoint of the walk, are functions of the adsorption activity a. In the desorbed 
phase (for a < a+) it is expected that R 2 ~ n 2 ", and H n ~ n u , where u — | is the metric 
exponent |14|. This scaling changes in the adsorbed phase (when a > a+) to R 2 ~ n 2 
and H n ~ constant. These expectations are confirmed by the data, as seen, for example, 


in figure [12] where data for the mean square radius of gyration are normalised and then 
plotted as a function of a. These graphs clearly show two scaling regimes, namely a 
high temperature phase (when a < a+) where the walk has bulk critical exponents and is 
desorbed, and a low temperature phase where the walk stays near the adsorbing boundary 
and has critical exponents of a linear object. 
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Figure 12: The mean square radius of gyration as a function of a, for adsorbing walks in Lj. In the left 
panel R 2 is divided by n 2v . For a < a+ this shows that R 2 ~ n 2u , and for a > 3+ R 2 increases in size 
faster than n 2 " (since it is adsorbed, it stays in the vicinity of the adsorbing boundary, and R 2 ~ n 2 in 
this regime. This is seen in the right panel, where R 2 is divided by n 2 . For a > at the curves collapse to 
a single underlying function, exposing the scaling in the adsorbed phase. The lengths of the walks varied 
from n = 50 (yellow) with colours increasing in hue to black when n = 500 in steps of 50. 


In general, the metric exponent associated with R 2 is a function of a, and it will be 
denoted by u a , where u a — | in the desorbed phase, and u a — 1 in the adsorbed phase. 
This exponent may be estimated from the mean square radius of gyration R 2 n data by 
examining the ratio 

log (Rl/Rl) 




log (n/m) 


(62) 


Here, v n m (a) is a function of n and m. By averaging over m, the estimate iv n (a) — 
(u n ,m(a)) m may be determined. In particular, fixing n and taking the average over m 
for 100 < m < 500 in multiples of 5 (and for m ^ n) gives estimates of iv n (a). These 
results are plotted in figure 13 for n £ (50,100,150,..., 500}. The data for a < 1.5 
gives v ~ 0.747, and for a > 1.95, u ~ 1.01. These results are evidence for the exact 
value v a — | in the desorbed phase, and v a — 1 in the adsorbed phase. 

The function u n (a) should scale with the combined variable r = n 4> (a — at). That 
is, one may expect that iv n (a ) = /7(r), where Id is a scaling function. In figure [14] the 
data in figure [13] are rescaled by plotting against r to uncover the scaling function Id. 

The (normalised) average height of the endpoint of the walk is plotted as a function 
of a in figure [15} These data show a clear transition where the scaling of H n changes. 
That is, the metric exponent associated with H n is a function of a, and is denoted by 


uf , where uf = v — | if a < at 


and ut — 0 if a > ai. 


so that H n ~ n u 


The graph 


of n 3 ^H n contains a set of curves which decreases with increasing a. These curves 
intersect each other close to at. 


An approach similar to the ratio method in equation (62) may be used to estimate 
iv (which will be referred to as the vertical metric exponent ): 

log (H n /H m ) 




g(a) = 


log(n/m) 


(63) 
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Figure 13: The metric exponent as a function of a, for adsorbing walks in L.f. The exponent was estimated 
from equation (62) by fixing n and then averaging over m for 100 < m < 500. With increasing n the 
change from v m | in the desorbed phase, to v = 1 in the adsorbed phase will become a step function 
at a = a+. The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black 
when n = 500 in steps of 50. 


This estimate of vj- is a function of n and m, and may be averaged over m to obtain 
iT-(a) = (u nm {a)) n . Taking the average for 100 < m < 500 in multiples of 5 to 
estimate v^(a) gives the curves in figure 16 when plotted against the combined variable 
r = n^(a — a+). It is seen in the graph that if a < a+, then u n (a) ~ but for a > a+, 
H n ~ const so that u^-(a) ~ 0. 



Figure 14: The same data as in figure 13 but now plotted against the rescaled variable r = n^(a — a+), 


for adsorbing walks in Lj. This rescaling of the activity a collapses the curves to a single underlying curve 
characterising the metric scaling of the walk in the desorbed and adsorbed phases. The lengths of the 
walks varied from n = 50 (yellow) with colours increasing in hue to black when n = 500 in steps of 50. 
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0 12 3 a 5 6 


Figure 15: The mean height H n of the endpoint of adsorbing walks in Lj of length n as a function of a. 
On the left n~ u H n is plotted as function of a. This quantity decreases sharply as a passes through the 
adsorption critical point The lengths of the walks varied from n = 50 (yellow) with colours increasing 
in hue to black when n = 500 in steps of 50. 


Since ~ n 2i/a where u a — u — | if a < a+, and u a 
may be defined by 


1 if a > a+, ratios of R 2 


n ^ o-2i/ 3 

Rln^ 


or 


2 2u R 2 

~RL 


l, 

—2t/ a _ 2 -t/2 


if a < a+; 
if a > a+. 


(64) 


0 2 2ur2 

the quantity is plotted as a function of the rescaled variable r = 
n 1/2 (a— a+) for n from 25 to 250 in steps of 25. In the desorbed phase this ratio should 



Figure 16: The vertical metric exponent v 1 as a function of a for adsorbing walks in Lj. The exponent 
was estimated from equation (63) by fixing n and then averaging over m for 100 < m < 500 . With 
increasing n the exponent changes from = v = | in the desorbed phase, to W = 0 in the adsorbed 

phase. Data for n = 500 and t > 5 did not converge, and is left out of this figure. The lengths of the 
walks varied from n = 50 (yellow) with colours increasing in hue to black when n = 500 in steps of 50. 
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Figure 17: Metric scaling for adsorbing walks in Lj. (a) The ratio R^n/{n 2l2 of the mean square radius of 
gyration plotted as a function of r = n^{a — a+) is displayed on the left. The ratio is given in equation 
([64]) and the data shows that it is approximately equal to 1 in the desorbed phase and approximately equal 
to 2 3 / 2-2 « 0.71 in the adsorbed phase, (b) The ratio H 2n /{n u H n ) plotted as a function of t. This ratio 
is defined in equation (65) and is approximately equal to 1 in the desorbed phase, and equal to 2 3 / 4 in the 
adsorbed phase. The data in both panels are for n e {25, 50, 75, 100,..., 250} and increase in hue from 
yellow (n = 25) to black (n = 250). 


be equal to 1, but, in the adsorbed phase, it should be equal to 2 1 5 0 0.71. This 

is clearly seen in the graph. The curves coincide well over the entire range of n and r, 
and decreases sharply close to the critical adsorption point at r = 0. A similar approach 
using the heights of the endpoint involves the ratios 


Hn 
H 


or 


2 n 


2 "Hn 

H 2n 



= 2 3/4 , 


if a < ; 

if a > a+. 


(65) 


The quantity 


2 "Hr 
H2n 


is plotted against r on the right panel in figure 


17 


In the desorbed 
n adsorbed phase the scaling 


phase this ratio should be equal to 1, as seen in the graph, 
of H n changes, and the ratio should be equal to 2 3//4 & 1.68, as seen in the graph. 


3.2. The generating function 

The generating function is given by the series 

oo n 

G(a, t) = ^^c+(v)a v t n . (66) 

n= 0 v=0 

Approximations of G(a, t) is given by the truncated sum 

N n 

G N (a,t) = J2J2 c "( v '>a v t n , (67) 

n=0 v=0 

where c+(v) is approximately enumerated by the GAS algorithm. In this study, G 50 o(a, t) 
was estimated using the approximate values of c,|(v) obtained in the simulations. 
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Figure 18: The phase diagram of the generating function G(a, t). The radius of converge of G{a, t ) is 
denoted by t+(a), and it is the critical curve. This curve is divided into two parts (a transition to a low 
temperature phase marked by To and a transition to a high temperature phase marked by A). The To and 
A curves meet in the critical point ft (located at the critical adsorption transition a+). The critical curve 
separates a "finite phase" (where G(a, t) is convergent and dominated by walks of finite length) from an 
“infinite phase" (where G(a, t ) is divergent). The scaling of G(a, t) is described by setting up a coordinate 
system (g, cr) with origin in ft, where o = (a—a+) and g = T — t. The scaling of G(a, t ) on approach to 
the critical curve is given by equation (69), depending on whether o < 0, o = 0, or o > 0. 


The critical curve t+(a) of G(a, t ) is its radius of convergence as a function of a. 
By equation (J8J) and by equation (J6]) , 

log £+(a) = -.F(a) = - lim £ log Z„(a). ( 68 ) 


The critical curve is shown schematically in figure 18 G(a, t ) is singular when t — £+(a), 
and if t > £+(a), then G(a, t) is divergent. 


The critical curve is parametrized by the scaling fields (cr, g) as shown in figure [18 
The critical point when a = a+ is located at (a+, £+), where £+ = t+(a+), and it divides 
the critical curve t+(a) into two parts. The part marked by t 0 corresponds to a transition 
to desorbed walks, so that £+(a) = l//x 2 in this regime (which is a transition to a high 
temperature phase). For a > a+ the approach to £+(a) is to adsorbed walks, along the 
critical curve marked by A (which is a transition to a low temperature phase). 

The phase diagram in figure [18] may be described in terms of the coordinates 
g — (^ — t) and cr = (a —a+). The behaviour of G(a, t) along its singular points 
along the critical curve is described by 

along T 0 ; 


G(a, t) 


-7l 


-7s 


-7+ 


at the critical point a = a + - 
along A. 


(69) 


In two dimensions exact values are known for the exponents: The exponent 71 = [ 101 

is the entropic exponent of half-space walks, and 7 S = || 


is the surface exponent 
of adsorbing half-space walks at the critical point a+. The exponent 7 + is the entropic 
exponent of adsorbed walks, and is given by the entropic exponent of self-avoiding walks 
in one dimension lower. In one dimension, this is 7 + = 1. 

The horizontal plane is the (a, £)-plane, 


A plot of G w (a, t) is shown in figure 19 


and the critical curve in figure 18 is shown as a black curve with the critical point shown. 
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Figure 19: A plot of logG 5 oo(a, t ) as a function of (a, t) for adsorbing walks in Lj- The critical curve is 
the black curve along the surface, and the critical point separating the To and A curves (see figure 18) 
is denoted. For t < t+(a) the generating G(a, t ) is finite, and for t > t+(a) it is divergent. These two 
regimes are clearly visible. 


Below the critical curve G(a, t) is finite, and approximated well by G N (a, t ) (for large 
values of N not too close to the critical curve). Above the critical curve G(a, t) is 
divergent, while G N (a, t), which is a polynomial, is finite. 

The exponent 71 can be estimated by putting a = 1 so that G(l, t) ~ <?“ 71 where 
g — (t+(l) — t) and £+(l) = ^ (and where g 2 is the growth constant of the walks in 
two dimensions). Thus, estimate 'y 1 by noting that 


log G(l, t) 
log g 


-71 + 


Cl 

log 9 


T 


C 2 

log 2 g 


+ .... 


(70) 


Proceed by approximating G(l, t) by G/v(l, t) (with t < £^~(1) so that ^ 5 oo(l. t) is a 
good approximation of G(l, £)). A least squares fit of the ratio on the left to a quadratic 
in gives the estimate 

7i = 0.952.... (71) 


This result is very close to the exact value || = 0.953 .... Similar analysis for a > 1 and 
a < 3 + gives results slightly larger, since the critical point at influences the data in 
its vicinity for finite values of N. 
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Figure 20: Scaling of the generating function near the critical adsorption point for adsorbing walks in Lj. 
The data points above are collected for g e [0.01,0.05] and y/g < a < y/g. The surface exponent 
7 S = , its exact value. 


A similar analysis with a — gives the estimate 

7 S = 1.429..., (72) 

for the surface exponent at the critical adsorption point. This is close in value to the 
exact result || = 1.453.... 



Figure 21: This figure verifies the scaling of the partition function for adsorbing walks in Lj as proposed in 
equation (77). The curves are plots of n _29 / 64 Z„(a)/U-T against the combined variable |r| = n 1 / 2 ^ —a+| 
(with the critical point approximated by a+ = 1.78 and a < a+)- The curves are walks of lenghts 
n e (50,100,150,..., 500}. All the data collapse to expose the scaling function h in equation (77). The 
lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black when n = 500 in 
steps of 50. 
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The case that a > a+ may also be analysed. The adsorbed walk should have the 
statistics of the self-avoiding walk in one dimension, so that the entropic exponent is 
7 + = 1 (this is the value of the exponent 7 in one dimension). Putting a — 3.5 and 
plotting Zn /n ( 3.5) against n gives an estimate for the critical value £+(3.5). In this case 
the data quickly converges to £+(3.5) = 0.260... (to three decimal places). Assuming 
that £+(3.5) = 0.260 and choosing the scaling field g — (0.260— £) gives the model 


log G(3.5,t) 
log ^ 


- 7 + + 


Ci 


+ 


C 2 


log g log 2 g 


+ 


(73) 


similar to the above, but now with the exponent 7 +. Plotting the left hand side as 
a function of k — j^, and fitting the data to a quadratic in k, give the estimate 
7 + = 1.00.... There remains, however, some curvature in the model for small values 
of g (and of cr), so that there may remain strong systematic corrections to this result 
(however, 7 + = 1 is consistent with this result). 

In the vicinity of the critical point (a+, £+) the generating function should exhibit 
scaling given by 


G(a, £) - g 7s f(g 0 cj), 


(74) 


where p — \ is the crossover exponent and f is a scaling function. That is, plotting 
g ls G(a, £) against the combined variable g~ <t> a should expose the scaling function f. 


n figure 20 this is done by plotting log(g 7s G(a, £)) against g 1 / 2 cr 


for g — — — f g 

_ ^ M2 

[0.01, 0.05] and cr = a —a+ £ [—\g 1/2 , g 1/2 ]. 

The partition function (see equation ([ 6 ])) should also exhibit scaling for large n, 
given by 


(75) 


where log/x a = £F(a) and log/x a = — log£+(a). The exponent 7 t can be related to 
the 7 -exponents as a —> a+ (that is, as cr —>• 0 ), and namely to the surface exponent 
7 S , where the walk is critical with respect to the adsorption transition. By noting that 
G(a, £) = ^ n Z„(a) t n , and approximating the summation by an integral, it follows that 
G(a, £) ~ g~ lt h( 0) if a = a+ and g — £+— £. This shows that 

7 t -l=7s-l = i-l = i. (76) 

That is, when a is close to a+ (and a < a+), the partition function has asymptotic 
behaviour 


Z n (a)^a 29 / 64 /7(^(a+-a)) M ". 

This result may be tested by plotting n' 29 / H Z„(a) (£+(a)) n against |t| = 


(77) 


- 1 -/ 2 I 7) — 3+ I 


All the data should collapse to the same universal curve exposing the scaling function h. 

3 + (and 


This is done with a < aj 


t+ 


and 0 < a < a+. 


[a) — Z;) in figure 


21 


for n £ {50, 100, 150_ 500} 


3.3. Adsorbing walks in the cubic lattice 

The (finite size) free energy JF n (a) is a function of the combined variable r = n 1 / 2 (a — 3 +) 
(see equations ( 10 ) and (fl 2 j); note that <fi — | for adsorbing walks in three dimensions 
1191). Plotting fF n(a) against r for data in the cubic lattice gives a graph similar to figure 
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Figure 22: The scaled free energy n{T n {a) — log fid) as a function of r = n 1 / 2 (a — a+) for adsorbing 
walks in Lj. The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black 
when n = 500 in steps of 50. The data collapse for small |t| (this is the critical scaling regime), and will 
approach a limiting curve as n —> oo, also for large r, given by the limiting free energy. 


[6} In figure 22 the scaled free energy n(fF n (a)— logyx 3 ) is plotted against r. The data 
shows a clear transition in the model from a desorbed to an adsorbed phase. 

Derivatives of the free energy to log a gives the (finite size) energy density £ n (a) 
and (finite size) specific heat C n (a). These are plotted in figure 23 against log a and in 

the critical point was approximated 
by = 1.31. This is a close approximation of the best estimate for the critical point 


figure [24] against r. In these plots, as in figure) 
by = 1.31. This is a close appn 
from our data (see equation ([79])). 



-3 -2-10 1 | 0 g 3 3 



i-1-1— \ - 1 — ~ — 

-3 -2-10 1 | 0 g a 3 


Figure 23: The energy £„{a) and specific heat C n (a) (see equation (fl3l)) as a function of log a for adsorbing 
walks in the cubic lattice. See figure [7] for the similar results in the square lattice. The lengths of the 
walks varied from n = 50 (yellow) with colours increasing in hue to black when n = 500 in steps of 50. 
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-30 0 30 r 60 -30 0 30 t 60 

Figure 24: The scaled energy density i/n£„(a) (left) and the specific heat C n {a ) (right) for adsorbing 
walks in Lj as a function of the r = n 1 i 2 {a — a+) for n in steps of 50 to 500. The lengths of the walks 
varied from n = 50 (yellow) with colours increasing in hue to black when n = 500 in steps of 50. The data 
collapse for small values of |r| (this is the critical scaling regime which contains the critical point). For 
larger values of |r| the curves will approach a limiting curve as n —> oo given by the limiting free energy. 


3.3.1. Location of the critical point a 


can be determined using the same analysis as in section 3.1.1 


The location of the critical adsorption point 

and in particular using 


equation (50) as a starting point. That, is for given values of n and m, an estimate 


+ 


Of 3+ 


can be obtained by solving for a in 
log(n£„(a)) log m 


log(m£ m (a)) log n ^ 

Here, the choice m — n— 100 worked well, and n was assigned values starting at n — 200 
to n — 500 in steps of 1. The estimates a+ = 3+„_ 100 showed a dependence on n, 
systematically decreasing with increasing n. The best estimate is obtained by using 
the model at — at — - 

n c | 0 

n G [200, 500] gives the best estimate 


suggested by equation (51). A least squares fit for al 


at = 1.3055 ± 0.0061. 


(79) 


The confidence interval is obtained by doubling the square root of the variance of the 


estimates a+. This result is slightly smaller than the result in reference 


namely 


= 1.334 ± 0.027 (obtained by using a Multiple Markov Chain implementation of the 
Berretti-Sokal algorithm |6|). 

The estimate (79) can be used to determine the crossover exponent 0. This is 
again done by considering the scaling of the specific heat (equation (JT3J)). It is expected 
that C n {at ) ~ n a(t, h c ( 0) when a — a+. An estimate of a<p is obtained by computing 
l09 |og^ c ~ f° r a ran 9 e °f values of n (in this case 100 < n < 500). The average is taken 
as the best estimate and a confidence interval is estimated by doubling the square root 
of the variance of the estimates. This gives oup — 0.0106 ± 0.0116. Determine the best 
estimate for 0 by using equation ( [TTj ): 

0 = 0.5053 ± 0.0053. 


(80) 
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Figure 25: The specific heat curves in figure 
critical adsorption point 


23 


for adsorbing walks in Lj intersect each other close to the 
This figure is a magnification of the curves close to where they intersect. See 


fig Lire [9] for the similar results in the square lattice. The lengths of the walks varied from n = 50 (yellow) 
with colours increasing in hue to black when n = 500 in steps of 50. 


This result compares well with the estimate 0 = 0.5005 ± 0.0036 for adsorbing walks in 


reference 1231. 


3.3.2. The critical point a+ and the specific heat C n (a ): The best estimate for a+ above 


(see equation (79)) should be examined by comparing it to estimates obtained from the 
specific heat curves in figure |23j These curves intersect each other near a+, and the 


region containing the intersections is magnified in figure 

The locations of the intersections between the specific heat curves in figure [25 
a function of n and are estimates of the critical adsorption point a\ 
intersections between C n (a ) and C n+W0 (a) against ^ (see figure 


is 


£ 


26 


By plotting the 
where n — 2 N 

and N e [23,200]); it is seen that the intersections fall approximately along a curve, 
which may be extrapolated to its intersection with the vertical axis. This gives a rough 
estimate of the critical point log a+ as being located in the the interval [1.33,1.35]. A 
more accurate extrapolation is done by using a linear least squares model to extrapolate 
to n — oo. Fitting to the model 3+ + ^ + ^, for all n > 50, gives the estimate 
log ~ 1.337. By examining the spread of the data in figure 26 a confidence interval 
can be estimated: 


log a+ = 1.337 ±0.020. 


(81) 


This estimate is slightly larger than the estimate in equation ( [79] ), but is consistent with 
the estimate a+ = 1.334 in reference |23|. However, the noise in the data in figure 26 
makes this a less reliable estimate. 


Equation (56) is equally valid for adsorbing walks in the cubic lattice. The last term 
on the right hand side is equal to zero when a — a+. Thus, by plotting log (C n (a)/C m (a)) 
against log(n/m), a set of curves should be seen which intersect when a — a+. At this 
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Figure 26: Extrapolating the intersections between specific heat curves for adsorbing walks in Lg . These 
data points are the locations of intersections between specific heat curves for n and for n+100 (and for 
even n and n e [46,400]). The points (except for some points at the largest values of n) approximately 
line up linearly if plotted against 1 fyfn. The location of the critical point can be estimated by extrapolating 
the line to the left vertical axis; this gives the estimate in equation ((TiT 


point the coefficient of log ^ is an estimate of acp. Since a — 0 in this model (and 
<fi — |), the critical point can also be determined by solving for a in 

'09 (£$) = °- < 82 > 

Solving this for n e [150,500] and m e [n— 100, n — 10] gives a large collection of 
estimates. The average is 

a+= 1.324 ±0.012, (83) 


where the confidence interval is one-half of the largest difference between two estimates 
in the collection. This result is smaller than the estimate in equation (81), and larger 
than the best estimate in equation (79). 

These results indicate that there may remain sources of systematic errors in the data 
and in the analysis, and that the estimates for a+ should be considered in this context. 


3.3.3. The microcanonical density function: The microcanonical density function of 
visits in adsorbing positive walks is determined from the microcanonical data in the 
model, and is given by 

P + (e)— lim (c+(|_e/ij)) 1/n = lim P+(e), (84) 

n —^oo n—>oo 

where P+(e) = (c+(LenJ)) lAl is a finite size approximation to P + (a). Existence of P + (e) 
can be shown (see for example reference |22|), and log P + (e) is a concave function of e. 

P + (e) can be determined by interpolating the finite size approximations P n + (e) and 
then extrapolating to n = oo by fitting a least squares model to the data. In figure [27] 
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Figure 27: The microcanonical density function P + (e) of visits for adsorbing walks in L 3 . These data 
the finite size approximations P+(e) to P + (e) for n = 100 and n = 500, as well as the extrapolated 
estimate of P + (e). The right derivative of log P + (e ) at e = 0 is an estimate of the location of the critical 
adsorption point a+ in the model. 


the data for the extrapolated function P + (e) is plotted together with P+(e) for n — 100 
and n = 500. 

A least squares fit of a quadratic to logP + (e) for e e [0,0.1] gives logP + (e) « 
1.54378 — 0.28704e —0.08784e 2 . By taking the right derivative and then taking e —» 0 + , 
an estimate for the critical point is obtained: 

a+ « 1.332. (85) 

The free energy is the Legendre transform of logP + (e). This may be estimated by 
fitting a polynomial to logP + (e). If a cubic polynomial in e is fitted to logP + (e) for 
0 < e < 0.5, then the estimated free energy for a > 4 is approximately 

F(a) « 1.6526 - 0.4240 log a - (0.3184-1.319 log a)y/~ 0.2438 + 1.0105 log a. (86) 

The critical point can be estimated as that location where the square root in the above is 
zero. This gives 3 + ~ 1.273. Similarly, the factor (0.3184 — 1.319 log 3 ) vanishes when 
at ~ 1.273. These estimates are far less secure than the estimates above, and are also 
smaller. 


3.3.4. Metric data: The mean square radius of gyration P 2 and mean height H n of the 
endpoint of the walk can be calculated as a function of a. In the desorbed phase (for 
3 < at) it is expected that P 2 ~ n 2u , and H n ~ rP, where u = 0.587... 111] is the 


metric exponent. This scaling changes in the adsorbed phase (when a > 3 +); in this 
phase it should be the case that P 2 ~ n 3 / 2 and H n ~ constant, since an adsorbed walk 
in the cubic lattice should have the statistics of a walk in one dimension lower. 


These expectations are supported by the data, as seen, for example, in figure 28 


where data for the mean square radius of gyration are plotted as a function of a. These 
graphs clearly show two scaling regimes, namely a high temperature phase (when 3 < 3 +) 
where the walk has bulk critical exponents and is desorbed, and a low temperature phase 
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Figure 28: The mean square radius of gyration as a function of a for adsorbing walks in L 3 . In the left 
panel R 2 is divided by n 2tJ . For a < a+ this shows that R 2 ~ n 2u , and for a > R 2 increases in size 
faster than n 2 " (since it is adsorbed, it stays in the vicinity of the adsorbing boundary, and it should be 
the case that R 2 ~ n 3 / 2 in this regime. This is seen in the right panel, where R 2 is divided by n 3 / 2 . For 
a> a+ the curves collapse to a single underlying, exposing the scaling in the adsorbed phase. The lengths 
of the walks varied from n = 50 (yellow) with colours increasing in hue to black when n = 500 in steps of 
50. 

where the walk stays near the adsorbing boundary and has critical exponents of a walk 
in one dimension lower. 



Figure 29: The metric exponent as a function of a for adsorbing walks in Lj. The exponent was estimated 
from equation (87) by fixing n and then averaging over m for 50 < m < 500. With increasing n the change 
from v = 0.588... in the desorbed phase, to u = 0.75 in the adsorbed phase will become a step function 
at a = a+. The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black 
when n = 500 in steps of 50. 
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Figure 30: The metric exponent as a function of the rescaled variable r = n^(a— a+) for adsorbing walks 
in Lj. The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black when 
n = 500 in steps of 50. 


The metric exponent u may be estimated from R* by examining the ratios 

2^(a) = '° g( f"; f? t ) - (87) 

log (n/m) 

Here, u nm (a) is a function of n and m. By averaging over m, the estimate iv n {a) — 
Wn,m) m can b e determined. Taking the average for 100 < m < 500 in multiples of 5 
(and for m ^ n ) gives an estimate for u n (a). The results are plotted in figure 29 for 

n G {50,100,150. 500}. The data for a < 1.2 give u ~ 0.592, and for a > 1.7, 

u ss 0.740. 

The scaling of iv n (a) as a function of r = n^(a —a+) can be uncovered by plotting 
the data in figure [29} This gives a set of curves which are very close to one another, 
uncovering a scaling function Id where u n {a ) = R{r). 

The average height of the endpoint of the walk is plotted as a function of a in figure 


31] The left panel displays the height normalised by division with n v and gives a set of 
curves which increase with n for a < a+, and decrease with n for a > a+. The curves 
intersect close to the critical adsorption point, and the limiting curve (in the n —> oo 
limit) should be a step function with critical point at a = a+. 

The vertical metric exponent u 1 - can be estimated from H n , by using a method 
similar to equation (87), namely an approximation by examing the ratios of /7„: 

log (H n /H m ) 




&) = 


( 88 ) 


log W m ) 

The exponent is approximated by n^(a) — (/x„ im (a)) n . Taking the average for 100 < 
m < 500 for fixed n gives esimates for u^-(a). If a < a+, then n^(a) should have value 
approximately equal to u\ that is, v n {a) ^ 0.58..., and for a > a+, H n ~ const so that 
//{-(a) S3 0 in this phase. The results are plotted in figure 32 against r = n^(a—a+). 
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0 12 3 p 5 6 


Figure 31: The mean height H n of the endpoint of adsorbing walks of length n as a function of a for 
adsorbing walks in Lg . On the left n~ u H n is plotted as function of a. This quantity decreases sharply as 
a passes through the adsorption critical point a+. This behaviour is also seen on the right, where n _1 /-/ n 
is plotted as a function n. The lengths of the walks varied from n = 50 (yellow) with colours increasing 
in hue to black when n = 500 in steps of 50. 


Since ~ n 2Ua where u a = u ^ 0.588... if a < a+, and v a — | if a > a+, the ratio 
of R\ n and R 2 n is given by 


R2 

n 2 _2t/ s 


Rin 


or 


2 2 "R 2 

R2 
' 'n 


1, 

2 2 ^—^ 2 ~ 0.32 


if a < 
if a > a+. 


(89) 



Figure 32: The vertical metric exponent as a function of a for adsorbing walks in Lj. The exponent 
was estimated from equation (88) by fixing n and then averaging over m for 100 < m < 500. With 
increasing n the exponent changes from W = u « 0.59 in the desorbed phase, to W = 0 in the adsorbed 
phase. The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black when 
n = 500 in steps of 50. 
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n" n 
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Figure 33: The rescaled mean square radius of gyration R 2 / n 2 " a for adsorbing walks in Lg plotted as a 
function of t = n^(a — a+) is displayed on the left. The metric exponent u a is the step-function defined 
by v a = f if a < a+ and v a = 1 if a > a+. The rescaled curves collapse to a single underlying universal 
scaling function. On the right the similar approach is used to uncover the scaling function for the mean 
height of the endpoint of the walk, scaled by The data are curves for n e {50,100,150,.... 500} 
increasing in hue from yellow (n = 50) to black (n = 500). 


n figure 33 this is plotted as a function of the rescaled variable r = n^ 2 (a — a+) for n 


from 25 to 250 in steps of 25. The curves coincide well with increasing n and signals 
a transition when a = a+ from the desorbed scaling regime into the adsorbed scaling 
regime. When a > a+ 2-°- 32 = 0 .801..., as shown in the graph. A similar approach 
using the heights of the endpoint would involve plotting 


H r . 




2" Hr 


1 , 


H 


or 


2/1 


H 


2 n 


->0.588 


if a < a+; 
if a > a+, 


(90) 


where iW is the vertical metric exponent. When a > a+ 2 0 - 588 = 1.503..., as shown in 


figure 33 


3.3.5. The generating function: The generating function of adsorbing walks in the cubic 


lattice is given by equation (66), where c+(v) is again the number of walks from the origin 
of length n in L 3 , and with v visits to the adsorbing boundary 9L 3 . Approximations to 
G(a, t) are given by G w (a, t) in equation (67), and G 50 o(a, t ) was calculated using the 
approximate values of c+(v) obtained by sampling with the GAS algorithm. The critical 
curve is given by equation (68) (see figure 18). G(a, t) is singular when t — t+(a), 
and if t > t+(a), then G(a, t) is divergent. In figure 34 the approximation G 50 o(a, t) 


is plotted, with the location of the critical curve, and critical point (a+,t+) indicated 
(where, as before, t+ = t+(a+)). The critical curve is similar to the critical curve in 
figure[l8j and the critical point divides the critical curve into two curves, namely a curve 
where the transition is a high temperature or desorbed walk marked by r 0 , and a curve 
where the transition is to a low temperature or adsorbed walk marked by A. Along r 0 
the critical curve is given by t+(a) = for a < and where g 3 is the growth constant 
of self-avoiding walks in the cubic lattice. 
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Figure 34: G 5 oo(a, t ) as a function of (a, t) for adsorbing walks in L 3 . The critical curve t c (a ) is the 
black curve, and the location of the critical point at a = a+ is indicated. Below the critical curve G(a, t ) 
is finite, and approximated well by G/v(a, t) (for large N and not too close to the critical curve). Above 
the critical curve G(a, t) is divergent, while G w (a, t ) is finite. 


The critical curve is parametrized by the scaling fields ( o,g ) as shown in figure 18 
Here, the scaling fields are given by g — t+ —t, and cr = a — a+. The singular points in 
G(a, t ) along t+(a) are described by the scaling assumptions shown in equation (69). 
The exponent 71 can be estimated by putting a = 1 so that G(l, t ) ~ g“ 71 . By using 
the model in equation (70), the estimate 

7 i = 0.725... (91) 

is obtained. This result is close to the estimate 0.697(2) in reference 1191. A similar 
analysis with a = a+ gives the estimate 

7 S = 1.203... (92) 

for the surface exponent at the critical adsorption point. This is slightly smaller than the 
estimate 1.304(16) in reference |31|. 

The situation is less clear in the adsorbed phase. The adsorbed walk should have 
the scaling of a self-avoiding walk in two dimensions, so that 7+ is given by the entropic 

exponent of walks in d — 2 (7+ = |§). Putting a = 4 and plotting Z}/ n ( 4) against n 

gives an estimate for the critical value t+(4). In this case the data quickly converges 
to t+(4) = 0.0926... (to four decimal places). Assuming that t+(4) = 0.0926 and 
choosing the scaling field g — (0.0926— t) gives a model similar to equation ([73]). 
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Figure 35: Scaling of the generating function near the critical adsorption point for adsorbing walks in L}. 
The data points are above were collected for g e [0.01, 0.05] and \yfg < a < \y/g. The value of the 
exponent g s was assumed to be given by the estimate in equation (76). 


Plotting the left hand side as a function of and fitting it to a quadratic for 
t e [0,0.06], gives the estimate 7 + ~ 1 . 2 ..., still well below the expected result 
|| = 1.34375. Examination of the data shows strong dependence of this result on 
the range of g in the model. For example, a fit with t e [0, 0.09] gives a smaller value 
7 + « 1.1.... These variable results indicate that G 500 (a, t ) is not a good approximation 


to G(a, t) near the critical curve A in figure [18] for adsorbing walks in the cubic lattice. 

In the vicinity of the critical point (a+, £+) the generating function should exhibit 
scaling given by equation ( [74] ). Plotting g ls G(a,t ) (with 7 S = 1.203) against the 
combined variable g~ 1 ^ 2 o should expose the scaling function f. In figure 35 this is 


done by plotting log(g 7 s G(a, £)) against g 1/2 o for g — — t e [0.01,0.05] and 

ct = (a-a+) e [ -\g x/2 , \g 1,2 \- 

The partition function (see equation ([ 6 ])) also exhibit scaling for large n. The 
scaling assumption in equation (75) applies here as well, where \ogg a — JF(a) and 
log /U a = — log £ c (a). As before, the exponent 7 1 is related to the 7 s -exponent as in 


equation (76): 


7 t - 1 = 7s - 1 « 1.203 - 1 = 0.203. 


(93) 


That is, when a is close to a+ (and a < a+), then the partition function has asymptotic 
behaviour 


Z„(a)~n°- 2m h(n*(at-a))n1. 


(94) 


This result may be tested by plotting a°- 203 Z n (a+) (t c (a)) n against r = a^(a — a+). This 


scaling is seen in figure 36 for a < ^ and n e {50,100,150,... 500}; all the data 


accumulate along a single curve, exposing the scaling function h. 
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Figure 36: The scaling of the partition function as suggested by equation (94) for adsorbing walks in Lj. 
The curves are plots of n _0 ■ 203 Z n (a)/x - ' , as a function of |r| = n x l 2 \a — a+| (with a+ approximated by 
1.33 and a < a+). All the data collapse to a single curve, which is the scaling function h in equation 
(94). The lengths of the walks varied from n = 50 (yellow) with colours increasing in hue to black when 
n = 500 in steps of 50. 


4. Conclusions 

The adsorbing self-avoiding walk is a classical model in rigorous and numerical statistical 
mechanics, and have received considerable attention in the physics and mathematics 
literature 112, 

In this paper the feasibility of collecting data in the microcanonical ensemble on 
adsorbing walks using a flat histogram implementation of the GAS algorithm 1241 was 
considered. This is an approximate enumeration algorithm, and the data can be used 
to directly estimate partition and generating functions, from which thermodynamic 
functions such as the free energy and specific heat can be determined (see, for example, 
reference |22j). 

The implemementation of the algorithm was done using endpoint elementary moves 
on half space self-avoiding walks, and the algorithm sampled from a flat histogram with 
reasonable success in both length and energy in the square and cubic lattices. Analysis of 
the data gives good results, better than previous Monte Carlo simulations in, for example, 
references 119, 23], but not as good as exact enumeration data in references |3,4|. A 
significant advance of this algorithm is that its produces a large amount of microcanonical 
data. Modifications to obtained data with respect to other quantities are trivial, and can 
easily be implemented. The simulations reported here were done on a Dell Inspiron 
530 desktop machine, but note that the algorithm can be implemented in parallel on a 
cluster with each cluster generating an independent sequence. This should give radically 
improved statistical data. 

The success of the implementation suggests that this numerical method may be used 
on other models (collapsing self-avoiding walks |30,39], for example). However, it may 
be necessary to extend the method by introducing, in addition to the sets of parameters 


18 19 23 
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denoted by {Pe, u } and { 7 ^}, additional sets of parameters which are conjugate to classes 
of elementary moves. For example, in the model of collapsing walks (see figure [4]), 
the energy may be changed by A u e {— 2d+l, —2d+2 ,..., 2d — 1 } by an elementary 
move, and parameters may be introduced for each value of A u to achieve flat histogram 
sampling (in a way similar to the introduction of 7 ^ for elementary moves increasing the 
energy of the walk. This will increase the complexity of the implementation, but with 
the result that flat histogram sampling will be easier to achieve. 

In the models of square and cubic lattice adsorbing walks, the algorithm produced 
data which gave good estimates of the locations of the critical adsorption point. The 


best estimates are obtained from equations (53) and (79), namely 


4 = 


1.779 ± 0.003, 
1.306 ± 0.007, 


in the square lattice; 
in the cubic lattice. 


(95) 


These results can be used to estimate the crossover exponent <p association with the 


adsorption transition, and our best esimates are seen in equations ( |54| ) and ( |80| ): 

0.496 ± 0.009, in two dimensions; 

0.505 ± 0.006, in three dimensions. 


(p = 


(96) 


In addition, other quantities from which 3 + and <p can be estimated were examined, and 
results largely consistent with the above values were obtained (see, for example, equations 


(55), (58) and (60) for square lattice results, and equations (81), (83) and (85) for cubic 


lattice results). These numerical estimates are in good agreement with those presented 
in reference |23|, and also in reference |3| in the case of the square lattice. The estimate 
in this reference, obtained from exact series data, namely a+ = 1.77564, agrees with 
the estimate above to two decimal places. In the cubic lattice the estimate for 3 + above 
is slightly smaller than the estimates 3+ = 1.338 ± 0.005 |31| and 3+ = 1.334 ± 0.027 
in reference |23| (rounding up of this last error bar gives a confidence interval which 
includes 1.306). 

The signature of the adsorption transition in the metric quantities of the model was 
also examined. The scaling of these quantities with r = n^(a —a+) were plotted in figures 


14 16 and 17 in the square lattice, and in figures BOl 32 and 33 in the cubic lattice. 


These results show a transition strongly characterised by changes in metric scaling and 
verify the value of the metric exponent and its finite size scaling through the critical 
point. 

Finally, the scaling of the generating function and partition function in these models 
were examined. Our results strongly supports the conventional properties of the model, 
and the values of the exponents { 7 i, 7 s , 7 +} estimated here are consistent with exact 
values and other estimates in the literature. 

The results in the square lattice are consistent with the exact values of 71 and the 
surface exponent 7 S , and the generating and partition partition function exhibit scaling 
consistent with the value of 7 S , as shown in figures[ 20 ]and[ 2 Tj In the cubic lattice our data 
gave the estimates 71 sa 0.725 and 7 S sa 1.203. These values are in addition to estimates 


elsewhere in the literature (see references 191, |31|), and although the estimates here 
may be improvements on previous esimates, they remain uncertain. However, scaling of 
























Microcanonical Simulations of Adsorbing Self-Avoiding Walks 


45 


the generating function in figure (35), and of the partition function in figure 36 


is some 


evidence that the esimate for 7 S is at least consistent with the scaling in the model. 
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